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Abstract. Quark-hadron (QH) duality in lepton scattering off nucleons is studied with the resonance quark model. It is shown 
that in the case of neutrino scattering off an isoscalar target the duality is simultaneously observed for charged and neutral 
currents xF^ , i*o , and xF? N weak structure functions. 

We demonstrate that the QH duality can be useful property for modeling structure functions in the so-called resonance 
region. As an example it is shown that combining relativistic quark model predictions with duality arguments allows a 
construction of the inclusive resonance F^ structure function. 

Keywords: Bloom-Gilman duality, nucleon structure functions, weak single pion production 
PACS: 25.30.Pt, 13.15.+g 



INTRODUCTION 

Historically quark-hadron duality was discovered as a re- 
markable relationship between hadronic scattering am- 
plitudes in s (low energy scale interaction) and t (high 
energy scale interaction) channels. 

The realization of the QH duality in inelastic electron- 
proton scattering was discovered by Bloom and Gilman 
|Q] (for review see J2j). It was observed that the ep 
inelastic VW2 resonance data scales as the deep inelastic 
scattering (DIS) one. In practise, it was shown that the 
resonance vW^ = F 2 data, considered as a function of 
ft)' = 1+W 2 /Q 2 scaling variable, is averaged by the DIS 
scaling curve. 

The Bloom-Gilman (BG) duality constitutes a kind of 
equivalence between the resonance and DIS structure 
functions. As was obtained by De Rujula et al. the 
duality phenomenon can be explained by doing the twist 
expansion of the F% p structure function. In low-g 2 , the 
leading twist of the F 2 ep dominates, while in high-g 2 
higher twist correction become more relevant. Thus the 
appearance of the BG duality means the suppression 
of the contribution from higher twists. De Rujula et 
al. showed also that the Nachtman variable ^(x,Q 2 ) — 

2x/ ^1 + 4-x 2 M 2 1 Q 2 ^j is a correct scaling variable 

to discuss QCD scaling violations, and consequently the 
proper scaling variable to study the BG duality. 

From the practical point of view the duality means 
equality between two integrals: 

/ d$F™(t;,Q 2 es )K d$F DIS (l;,Ql IS ), (1) 

where F res and F DIS denote the structure functions at 
resonance and DIS regions respectively. The ^ min , B, max 



are defined by the limits of resonance region, namely: 

min j max ^ ( Q-res ) ^max/min ) ■ (2) 

For the lower bound of the resonance region we take 
Wmin = M + m K , while W max will vary from 1.6 to 
2.2 GeV. The typical four-momentum transfer in the res- 
onance region for ep scattering is < Q 2 es < 5 GeV 2 . 
For inclusive vN scattering we consider < Q 2 es 3 GeV 2 . 
The characteristic four-momentum transfer for DIS is 
Qdis ~ 10 GeV2 - With such kinematical settings F D!S 
is integrated over a different kinematical domain, in in- 
variant, mass than F res . 

Recently the BG duality has been carefully studied by 
Niculescu et al. J4). The new precision resonance JLab 
data has been discussed together with the old measure- 
ments. The BG duality was quantitatively studied in a 
wide kinematical region: 0.3 < Q 2 res < 5 GeV 2 , 1.1 < 
W 2 < 4 GeV 2 . The results of this analysis confirmed the 
appearance of the BG duality but with accuracy depend- 
ing on the prominent resonance regions. 

The duality is also observed in electron-deuteron scat- 
tering [5 1, which is also evidence for the duality between 
electromagnetic neutron structure functions. 

It seems that analogously as for the eN scattering, the 
duality should be observed also for weak nucleon struc- 
ture functions, which are used to describe the inclusive 
vN interaction. In contrast to the eN scattering the ex- 
perimental evidence for the duality in vN scattering does 
not exist. The neutrino scattering data, collected in the 
resonance region, is still too imprecise to order to extract 
directly the resonance structure functions. On the other 
hand the knowledge of the DIS structure functions for 
inclusive vN scattering is also limited (see |6)). 

Given the lack of appropriate neutrino scattering data, 
the duality in vN scattering has been studied only with 
the phenomenological descriptions 018] [9), which have 
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FIGURE 1. In the left panel, plots of p structure functions are presented. The quark model predictions (solid line) are compared 
with the inclusive CLAS data 1 18 1 (full squares and triangles and open circles) and DIS scaling curve (dashed line). In the right 
panel the plots of the ratios |4](, computed for W max =1.6, 1.8, 2.0 and 2.0 GeV, are shown. 



been fine tuned to the old neutrino-deuteron bubble 
chamber scattering data. 

Assuming that the duality is a property of the electro- 
magnetic and weak nucleon structure functions can lead 
to the additional constraint on: the parameterization of 
the nucleon and resonance form factor^] as well as the 
nucleon inclusive structure functions. 

This talk presents an extension of the analysis pub- 
lished in Ref. (8). The theoretical framework is formu- 
lated in Section 2. In Section 3 we demonstrate that com- 
bining quark model predictions and duality arguments 
leads to the reasonable parametrization of the inclusive 
Ef p structure function in the resonance region. The sec- 
ond aim of the talk is to investigate the BG duality in 
neutrino scattering off nucleons. In the last section we 
shortly discuss the duality in lepton scattering off nuclei. 



THEORETICAL FRAMEWORK 

In order to compute the resonance structure functions we 
apply an updated 02 Rein-Sehgal (RS) model 02). The 
RS model lfl2l was devoted to describe the 1 % production 
induced by neutrino-nucleon interaction in a wide kine- 
matical region (up to W = 2.0 GeV). The description is 
based on the relativistic harmonic oscillator quark model 
03). Even thought the RS model is relatively simple and 
old it is still the reference description for the experimen- 



tal data analysis in modem experiments like K2K and 
MiniBooNE. 

In the RS model the vector and axial hadronic currents 
are described with the nucleon vector and axial elastic 
form factors but, as it was shown in Ref. ifTTl . apply- 
ing these form factors leads to the underestimation of 
the vector part and overestimation of the axial contribu- 
tion. In this analysis we use the vector and axial form 
factors proposed in Ref. [11]. The updated vector con- 
tribution agrees with the recent fits of the resonance-like 
electroproduction data (in the ^33(1232) resonance re- 
gion), while the new axial form factor was obtained as 
a simultaneous fit to the ANL and BNL In production 
data^] The updated RS model works efficiently and fits 
very well to the recent MiniBooNE measurements 04). 

The mechanism for the duality appearance in lepton 
scattering off nucleons is partially explained by quark 
model arguments 03). The RS model as an example 
of a quark model with strong experimental applications 
seems to be an interesting framework for phenomenolog- 
ical studies of the duality in vN scattering. 

Close & Isgur claimed that one of the necessary 
conditions for the realization of the duality is taking into 
consideration at least one complete set of resonances 
of each symmetry-type. Thus, for completeness of our 
analysis we consider the following list of resonances^] 
P 3 3(1232), Pn(1440), Di 3 (1520), 5n(1535), 5 3 i(1620), 
Sn(1650), P 3 3(1600), Di 3 (1700), D, 5 (1675), 
F 15 (1680),D 33 (1700), P n (1710), P l3 (1720), P 31 (1750), 



1 In Ref. |3 1 local duality in elastic ep scattering is discussed. The 
applications can be found in Ref. [101. 



2 The axial form factor is given by Eq. 52 from the erratum of Ref. II 11 . 

3 The interference terms between resonances are also included in the 
calculus. 




FIGURE 2. In the left panel the predictions of the resonance structure function computed with Eq.[6](thin solid line) are shown. 
We plot also the scaling function (thick dashed line), valence scaling function (thin dashed line) and resonance fit, given by Eq. 9 
of Ref. |4|, (dotted line). In the right panel the plots of the ratios Eqs.[4| and|5| computed for 1.15<W<3.9 GeV 2 , are shown (solid 
and dashed lines respectively). The results are compared with the data from 141 . 



A3(1900),S 3 i(1900),F 3 i(1910),F 35 (1905),F 37 (1950), 
F 33 ( 1920) , F n ( 1 990) , F 15 (2000) , F 35 (2000) , G n (2 1 90) , 
Gi9(2250). However, one needs to remember that for 
larger values of the hadronic invariant mass the har- 
monic oscillator quark model predictions deviate from 
the experimental measurements. Therefore extending 
the number of resonances in the discussion can give only 
a hint on understanding the duality. 

In order to quantitatively study the BG duality we 
introduce the ratio: 

R i (fi . Qres > QdIs) = "J^ ) ( 3 ) 

/ di;f DIS (i;,Q 2 D IS ) 

where the limits of integrals are defined by while 
ff IS and f[ es are computed with GRV94 PDF'sjtm and 
the updated RS model respectively. 

F{ p STRUCTURE FUNCTION 

In order to test our theoretical model we compare the 
theoretical predictions for inclusive F 2 P structure func- 
tion with the resonance CLAS data [18]. As can be seen 
in Fig. [T] the quark model predictions are systematically 
below the data but it is not a surprise, because the non- 
resonant background (NoBG) contribution is not taken 
into account. On the other hand, the experimental data, 
and theoretical predictions seem to lie along the valence 
pep,val curve ( see pjg j^j, 'rhjs property of the resonance 
data has been already studied in Ref. |3). 



The quantitative verification of our approach is done 
by computing the ratio (see Fig.[T| 

Rf{Q 2 res ) = R (F' p ,Q 2 res , Q 2 DIS = 10 GeV 2 ) , (4) 

where we consider five different values of <!;„„„, namely 
Wmax =1.6, 1.8, 2 and 2.2 GeV. It is seen that, above 
g 2 es = 0.5, the ratios Q computed for W max = 2 GeV 
and W max = 2.2 GeV slowly vary on Q 2 res (R e 2 p (Q 2 es > 
0.5) - 0.6 and R 2 p (Q 2 es > 0.5) - 0.7 respectively). The 
weak Q 2 es dependence of R 2 p (Q 2 es ) is evidence for the 
same scaling behavior of the resonance predictions and 
the DIS structure function, however, it is clear that some 
strength is missing to get a perfect duality. 

The scaling is broken below Ql es = 0.5 GeV 2 , where 
R e 2 p (Q 2 es ) rapidly tends to zercrj The argument which 
explains the mechanism of the duality violation (below 
Ores — 0-5 GeV 2 ) was given in iflTl : in the quark models 
the duality is strongly related to the dominance of the 
magnetic contribution, while at low Q 2 , the electric and 
magnetic multipoles have comparable strengths. 

In order to compare more qualitatively with the exper- 
imental measurements we compute the ratio (Eq. [4]i for 
1.1 < W 2 < 3.9 GeV 2 and compare it with the analogous 
one obtained by Niculescu et al. [4| (here the DIS con- 
tribution was predicted with MRS(G) LQ PDF's). Both 
ratios are plotted in right panel of Fig. [2] Our result is 
systematically shifted down with respect to the experi- 
mental data. 



4 Gauge invariance implies that as F^'^i^es ~ * ") — > 0; as a conse- 
quence, in the ratio of Eq.|4] the numerator decreases faster than the 
denominator. 



In Ref. (4) it was shown that if one assumes the 
BG duality for F^ structure function, then the reso- 
nance data can be approximated by the DIS-like F^ s pm 
parametrization. The obtained fit (on an average level) 
describes well resonance data. We use this parametriza- 
tion in our analysis. In Fig. [2] we plot the ratio 



d$F™- par ($) 



computed for 1.1 < W 2 < 3.9 GeV 2 . We see that the 
above function weakly depends on Q% s . Additionally, in 
almost all considered domains R'^ 1 '"' (Q 2 es ) is around 
0.6. It illustrates the fact that the quark model systemat- 
ically underestimates the experimental data, and it needs 
to be enriched by the nonresonant background contribu- 
tion. On the other hand this regular deviation between 
quark model predictions and the data allows us to pro- 
pose the effective structure function as: 



2 I'l/TU.-ti,^;!^ , (6) 

T}(€) = 2arctan((W -M-m^j 'GeV)/ 'tc, (7) 



where the K 



ep,val 



is the valence DIS contribution, while 



Tj(§) function is introduced to make a smooth transition 
from the threshold to the ^33(1232) resonance peak. It is 



clear the ratio 



' Qres i Qr 



10 GeV and the 



analogously computed ratio (Eq. |3j will tend to one. In 
Fig.[2|(left panel), we plot F^ p dual together with the CLAS 
data^he good agreement in the first and also second 
resonance region can be seen. 

The above approximation of the structure function 
seems to be naive in simplicity but very effective in the 
applications. In particular, it can be useful to apply to the 
Monte Carlo generators, dedicated to simulate neutrino- 
matter interactions. 



We consider the weak charged current (CC) and neu- 
tral current (NC) nucleon structure functions. The fol- 
lowing ratios are computed: 



RTiQtes) 

<(gL) 
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10 GeV 2 ) , (8) 



Qdis 1 



2^=10 GeV 2 



(9) 



R(xF^,Q 2 res ,Q 2 DIS : 



10 GeV 2 ) (10) 



(5) We found out that for W max = 1.8 GeV (see Fig. [3]) 
the above ratios, similarly as for ep scattering, weakly 
depend on Q 2 es . Once again it can be treated as evidence 
for appearance of the duality. 

It is interesting that the ratios of CC and NC structure 
functions are very similar. It is the result of the Si/ (3) x 
ST/ (2) and SU(6) symmetries of the structure functions. 
As a consequence if the duality is observed for the one 
channel, say CC, it must be also visible (on the same 
level) for the NC structure functions. This property can 
have applications in modeling of the cross sections for 
the neutral current neutrino-nucleus scattering. 

As was mentioned in first section, the predictions of 
the cross sections for neutrino-nucleon scattering are 
affected by the lack of knowledge of the axial structure of 
the nucleon. Therefore we present the ratios (|8fT0|i with 
1 a error coming from the axial form factor uncertainty 
(we apply the results of Ref. |[T9l . where the uncertainty 
of the axial contribution was discussed). It is seen (in 
Fig. p} that for the xF^ N and xFf N structure functions the 



BG duality is observed within la error, whereas for F~ N 
structure function the duality appears within 2a error. 

However, we remark that in the discussion of structure 
function uncertainties, one needs to remember that in our 
approach the description of NoBG contribution is quite 
approximate. It gives rise to a systematic error. In order 
to show the role of NoBG contribution we plot the ratios 
(Eqs.lSHTOll computed without NoBG. 



OUTLOOK 



NEUTRINO-NUCLEON INTERACTIONS 

It is well known that if the ^33(1232) resonance region 
is included in the discussion then the BG duality can not 
be simultaneously observed for the proton and neutron 
targets. It is the result of the isospin symmetry of the res- 
onance structure function and SU (6) symmetry property 
of the DIS structure function. Fortunately the BG duality 
(assuming some accuracy) appears for the weak structure 
functions of the isoscalar target . 

For the study of vN scattering we apply the updated 
RS model but the description is supplemented by includ- 
ing the nonresonant background contribution (we follow 
the effective description proposed in Ref. lfl"2"ll ). 



As was mentioned in the first section, the duality is 
observed in electron-deuteron scattering. But the scaling 
behavior of the resonance data for electron scattering off 
heavier nuclear targets, like carbon and iron |20|, is also 
observed. Here the duality between the resonance and 
the DIS nuclear structure function is visible even before 
averaging the data. 

In fact, the Fermi motion effect naturally "averages" 
the nucleon structure function over a wide kinematical 
region, but, there are also other nuclear medium effects 
like renormalization of the resonance properties, redistri- 
bution of the resonance prominent regions etc.. It makes 
the theoretical investigation of the duality in lepton scat- 
tering off nuclei very difficult and delicate. 




FIGURE 3. Plots of Eqs.[8][9] and[lO] computed with W max = 1.8 GeV are shown. The solid lines denote the full model predictions 
with NoBG. The computations without NoBG contribution are denoted by dashed lines. The thick and thin lines denote the ratios 
for CC and NC structure functions. The shadow areas denote the lo~ uncertainty. 



On the other hand, observing the duality for nuclear 
targets means that the nuclear corrections in the reso- 
nance and the DIS regions are comparable. But it has 
been not expected 0. 

Some effort to investigate the duality in lepton scat- 
tering off nuclei was done by Lalakulich et al. lETl . 
The electromagnetic and weak nuclear structure func- 
tions were obtained with an independent shell particle 
model, and then compared with the experimental data. 
For the v-nuclei scattering the duality was not obtained. 
However, probably further studies of the nuclear effects 
and the role of NoBG contribution can shed light on this 
topic [22]. 

It is obvious that assuming that the Bloom-Gilman 
duality is also a fundamental property of weak nuclear 
structure functions can give interesting applications in 
neutrino scattering physics. In particular in can be useful 
for: 

(i) modeling the neutrino cross sections in the reso- 
nance region; 

(ii) extraction of the axial contribution from the scatter- 
ing data; 

(iii) studying the nonresonant background contribution. 

Hopefully the forthcoming neutrino experiments, like 
Minerva will be able to critically investigate the duality 
in neutrino scattering off nuclei. 
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